Neutrino Constraints on the Dark Matter Total Annihilation Cross Section 
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In the indirect detection of dark matter through its annihilation products, the signals depend on 
the square of the dark matter density, making precise knowledge of the distribution of dark matter 
in the Universe critical for robust predictions. Many studies have focused on regions where the dark 
matter density is greatest, e.g., the Galactic Center, as well as on the cosmic signal arising from all 
halos in the Universe. We focus on the signal arising from the whole Milky Way halo; this is less 
sensitive to uncertainties in the dark matter distribution, and especially for flatter profiles, this halo 
signal is larger than the cosmic signal. We illustrate this by considering a dark matter model in which 
the principal annihilation products are neutrinos. Since neutrinos are the least detectable Standard 
Model particles, a limit on their flux conservatively bounds the dark matter total self-annihilation 
cross section from above. By using the Milky Way halo signal, we show that previous constraints 
using the cosmic signal can be improved on by 1-2 orders of magnitude; dedicated experimental 
analyses should be able to improve both by an additional 1-2 orders of magnitude. 

PACS numbers: 95.35.+d, 98.62.Gq, 98.70.Vc, 95.85.Ry 



I. INTRODUCTION 

Dark matter, an invisible substance originally pro- 
posed by Zwicky in the 1930s to explain the mass-to- 
light ratio of the Coma galaxy cluster [l[, still evades 
revealing its true identity after many decades. Subse- 
quent works have confirmed its pervasive existence, and 
its contribution to the critical energy density has been 
measured to be ft x ~ 0.3. The distribution of dark 
matter (DM) in the universe has also been thoroughly 
studied on a wide range of scales. Many DM candidates 
have been proposed, from axions to sterile neutrinos to 
weakly interacting massive particles (WIMPs); see re- 
cent reviews 0, H H|. WIMPs are especially popular 
since they are among the best theoretically motivated 
non-baryonic DM candidates and also capable of produc- 
ing the correct relic density. They can be detected indi- 
rectly through the astrophysical signatures of their an- 
nihilation products, detected directly in low background 
experiments underground, and produced at particle ac- 
celerators. 

For indirect detection, the DM distribution is essential 
for predicting the detectability of the signals, as the an- 
nihilation rate depends on the annihilation cross section 
and the square of the DM number density (the DM mass 
density divided bythe DM mass). The Galactic Center 
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intermediate mass black holes U9LI20L |2lj|22l|23i, dwarf 
satellite galaxies 0, [H, [26| , and nearby galaxies 12LI271]. 
i.e., places where the dark matter is strongly concen- 
trated, have been considered in the search for DM sig- 
nals. Additionally, the diffuse cosmic signal due to 
all distant concentrations of DM has been utilized as 
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In the Milky Way (MW), an attractive site for con- 
sideration is the GC, due to an expected increase in the 
DM density towards the centers of gravitational poten- 
tials [13, [3a. l39ll . and possibly a spike around the central 
black hole [l4fl5l[l6|. However, the inner gradient of the 
halo profile is the most uncertain part. Although simu- 
lations generally suggest a steep cusp towards the halo 
center, some observations, such as the high microlensing 
optical depth, support much flatter profiles [4{j. In fact, 
there is no direct observational evidence that any nearby 
galaxy has a cusped DM profile. The small spatial resolu- 
tion needed to settle the debate has not yet been reached 
in simulations. Another problem is astrophysical back- 
grounds. For example, astrophysical gamma-ray sources 
in the GC may limit the DM parameter space that can 
be probed with gamma rays [4l| . So it may be preferable 
to study other regions of the Milky Way, where both the 
signal and backgrounds are more certain. 

In this study, DM annihilation in the whole MW halo, 
as opposed to just the GC, is our primary focus. While 
the GC signal strongly depends on the uncertain cuspi- 
ness of the underlying halo profile, the signal averaged 
over larger regions of the halo is much less dependent on 
the details as long as the density normalization is chosen 
consistently with observed data, e.g., rotation curves 
This has been emphasized by e.g., Refs. [43l l44l [45l] . 

Using the annihilation signal arising from the halo 
means accepting a lower signal intensity, but the in- 
tensity of the astrophysical backgrounds is also lower. 
The halo signal is less subject to theoretical uncertain- 
ties compared to the GC or the cosmic signal, and can 
provide more robust results. Moreover, regardless of 
the parametrization, the halo always provides a minimal 
isotropic contribution that cannot be distinguished from 
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any truly diffuse cosmic signal. The halo isotropic com- 
ponent is increasingly dominant over the cosmic signal 
for flatter halo profiles. 

In indirect detection, regardless of the identity of the 
DM particle, one would expect the most stringent con- 
straint on the annihilation cross section from the most 
detectable annihilation products, and such gamma rays 
have long been sought. However, the branching ratios to 
produce gamma rays are small and uncertain, making it 
impossible to set a certain limit on the total annihilation 
cross section. Since neutrinos are much harder to de- 
tect, one may naively expect that they cannot yield any 
meaningful constraints. However, this is not the case. 

While both the gamma ray and neutrino flux predic- 
tions are very much dependent on WIMP model param- 
eters, one can avoid these complications by adopting a 
more conservative perspective, following Beacom, Bell 
and Mack (BBM) [46| . If we assume that only Standard 
Model (SM) final states are produced, e.g., purely sterile 
neutrinos are not considered, then neutrinos and antincu- 
trinos are the least detectable final states. Anything else 
will inevitably produce the much more visible gamma 
rays, leading to a stronger cross section limit. Thus, one 
can obtain a stringent upper limit on the total annihi- 
lation cross section by assuming that only neutrinos are 
produced as final states, and comparing the expected sig- 
nal to the well-measured and modeled flux of atmospheric 
neutrinos, which is the dominant background. 

We focus on the importance of the annihilation sig- 
nal arising from the halo in the context of obtaining 
model-independent bounds on the DM annihilation cross 
section, as described above. For very large annihilation 
cross sections, as in Ref. [13], halos would be significantly 
modified. BBM have shown that such large cross sections 
are not possible [46| , since then the cosmic signal, which 
depends on the time-integrated annihilation rate of all 
halos, cannot be accommodated in comparison to the at- 
mospheric neutrino spectrum. In this study, we instead 
use the MW halo signal, which is generally larger and 
more certain, as we show below, to set an even stronger 
limit on total annihilation cross section. 

We first summarize the halo profiles and the model for 
annihilating DM in the halo and in cosmic sources. We 
also identify the various components of the annihilation 
signal arising from the halo. To provide the most general 
results, we consider various profile shapes for the halo 
together with a dark matter model in which neutrinos 
are the only annihilation products. Next, we compare 
the neutrino spectrum arising from annihilations to the 
atmospheric neutrino spectrum, and arrive at our con- 
straints on the total annihilation cross section. Our study 
shows that even the full-sky halo signal would lead to over 
an order of magnitude improvement on the cosmic sig- 
nal limit on the cross section. If the broad enhancement 
towards the GC direction is taken into account, further 
improvement is possible while remaining secure with re- 
spect to the halo profile uncertainties. 
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FIG. 1: The dark matter density versus radius for the three 
profiles considered (Moore, NFW and Kravtsov in order of 
dotted, dashed and solid lines), based on Eq. The nor- 
malizations are chosen to match the large scale properties of 
the Milky Way. The vertical gray line indicates the solar circle 
distance, R sc ~ 8.5 kpc. See Table I for the other parameters. 



II. ANNIHILATING DARK MATTER MODEL 

Cold dark matter simulations show that DM assem- 
bles hierarchically in halos and successfully predicts the 
formation of large scale structure. The general form 
of the halo is expected to be independent of the total 
mass enclosed over many orders of magnitude. A useful 
parametrization that fits the results of the simulations is 
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where among the variables (a,/3, 7, r s ), 7 is the inner 
cusp index and r s the scale radius. In Fig. 1, we plot 
some commonly used profiles, the Navarro-Frenk- White 
(NFW) [H|, Moore and Kravtsov [H profiles, with 
the corresponding parameters listed in Table I. The nor- 
malizations are chosen so that the mass contained within 
the solar circle (R sc — 8.5 kpc) provides the appropriate 
DM contribution to the local rotational curves (see Q 
for details of a similar analysis). This yields the local 
DM density pnfw(Rsc) = 0.3 GeVciTr 3 for the NFW 
profile; for the others, see Table I. 

In general, the profiles agree on larger scales. However, 
the profiles in the inner regions are more uncertain, due 
to the spatial resolution limit of the simulations. Physi- 
cally, baryons typically dominate these regions, and their 
impact has not yet been fully taken into account in the 
simulations. The effect of baryons is complicated and can 
increase the inner cusp strength [H, [4i| as well as provide 
a mechanism for momentum transfer to DM particles, 
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which works to flatten cusps [50(. With these points in 
mind, we choose the halo profiles shown in Fig. 1 to illus- 
trate the impact of variations in the halo parametrization 
on the signal predictions. 

A. Annihilations in the Halo 

The calculation of the spectrum of DM annihilation 
products from the halo is well established (e.g., Ref. 0]). 
For annihilations, the intensity (number flux per solid 
angle) at an angle ip with respect to the GC direction is 
proportional to the line of sight integration of the DM 
density squared, or, 

1 ftma* 

J WO = „— a- / p 2 WR 2 sc ~2lR sc cosi> + P) M, 

ttscPsc JO 

which we choose to normalize at the solar circle (R sc = 
8.5 kpc) with p sc — 0.3 GeV cm -3 . Note that the pref- 
actor \/{Rscp1c) is an arbitrary scaling that is used to 
make J dimensionless regardless of the underlying DM 
density profile. The upper limit of the integration, 



(■max = y (R 2 MW ~ ^ V^sc) + R sc COS Ip , (3) 

depends on the adopted size of the halo, Rmw- However, 
contributions beyond the scale radius, typically about 
20-30 kpc, are negligible. It is also useful to define the 
average of J in a cone with half-angle ip around the GC 
that spans a field of view of Ail = 2n(l — cos ip): 

i r 1 

Jau = -r^ / J(ip ') 2tt d(cos V') • (4) 

A" Jcosj/> 

Then the average intensity of DM annihilation products 
from the field of view can be cast as 

^AS] _ (<?AV) „ RscPlc dN 

where dN/dE is the spectrum of annihilation products 
and m x is the assumed mass of the DM particle. The 
factor 1/2 accounts for DM being its own antiparticle and 
1 /47r is for isotropic emission. The thermal average of the 
annihilation rate is proportional to (<Jav), the product of 
the annihilation cross section and the relative velocity. 

In Fig. 2, we show the behavior of J(ip) (thin lines) 
and Jaq. (thick lines), as functions of ip and AQ/4tt, re- 
spectively, for the three representative halo profiles. Note 
that J diverges for very cuspy profiles. This cusp may 
be an artifact of simulations, which cannot resolve small 
spatial scales. To avoid such numerical divergences, we 
put in a flat core for all the profiles in at the innermost 
0.1°, which corresponds to ~ 0.015 kpc; our results are 
not sensitive to how the inner profile is regulated. 

One can expect very large enhancements for small an- 
gles ip around the GC. Many studies have therefore fo- 
cused on emission from the few degrees around the GC. 
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FIG. 2: Line of sight integration J(ip) as a function of the 
pointing angle ip (bottom axis) with respect to the Galac- 
tic Center direction for the three different profiles considered 
(Moore, NFW and Kravtsov in order of dotted, dashed and 
solid thin lines). Its average, Jasi, inside a cone with half 
angle ip around the GC as a function of the visible fraction 
of the whole sky, Af2/47r = (1 — cosi/>)/2 (top axis), is also 
presented (thick lines). Note that the left (right) side of the 
graph is presented in log (linear) scale in Af2/47r. 



However, while this takes advantage of the higher DM 
densities around the GC, it also strongly depends on the 
chosen profile. As shown in Fig. 2, the differences in the 
intensity due to different choices of profile can be orders 
of magnitude at small angles. We therefore focus on the 
signal arising within angles ip ~ 30° or larger around 
the GC. This is motivated for three reasons. First, the 
enhancement J~aq is much less sensitive to the chosen 
profile for larger ip. Second, the viewing angle ip to con- 
sider depends on the particle considered, and for neutri- 
nos, the angular resolution is limited by the kinematics of 
the detection reactions to be 66(E) ~ 30° x v / [GeV]/E. 
Thus, it is impractical to consider the arrival direction for 
neutrinos less energetic than several GeV. Third, choos- 
ing a large field of view ensures that there are adequate 
statistics for the detection of the signal and background. 

We now define three signals corresponding to different 
components of the halo: 

• Halo Angular: We consider a 30° half-angle cone 
for neutrinos more energetic than 10 GeV and a 
larger half- angle of ~ 3 x 86(E) for lower energies. 
A conservative value for this Halo Angular compo- 
nent inside ip ~ 30° is Jau 9 — 25 (below 10 GeV, 
when larger angular regions are used, this would be 
smaller) . 
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• Halo Average: When a large field of view Afl is 
considered, the average J^n is much less sensitive 
to the chosen profile, as shown by the differences in 
J avi in Fig.[5J We define the Halo Average, J Ave — 
5, as a canonical value for the whole sky (ip — 180°). 
Note that the differences due to different choices 
of profiles from this canonical value are small, less 
than a factor of 2. 

• Halo Isotropic: Another key point is that J{^p) 
never vanishes, even in the anti-GC direction (ip = 
180°), due to a minimal isotropic component of the 
halo signal that arises within our immediate neigh- 
borhood. This is present for all profiles and cannot 
be differentiated from the truly cosmic signal, ex- 
cept by its spectrum. We define this minimal com- 
ponent as Halo Isotropic, Ji so ~ 0.5, which yields 
the most robust constraints. 

The values for these three components are summarized 
in Table I. Finally, halos may have substructure which 
can contribute an additional boost to DM annihilation 
signals, e.g., Refs. [H, [5l|, [52j. To be conservative, we 
do not attempt to model any substructure, but note that 
it would only enhance all of the annihilation signals and 
hence improve the constraints that we present below. 

B. Cosmic Annihilations 

The calculation of the diffuse cosmic DM annihilation 
signal has been thoroughly discussed in the literature 
(e.g., Ref. 29]). The intensity of the cosmic diffuse signal 
is 

(oav) tf x pl c f dN(E') (1 + zffjz) 
dE 2 AirmlHoJ dE' h(z) ' 1 ' 

where p c — 5.3 x 10~ 6 GeVcm -3 is the critical den- 
sity, Hq = 70kms~ 1 Mpc -1 is the Hubble constant, and 
h(z) = [(1 + z) 3 fl x + Qa] 1/2 with n x = 0.3 (ignoring the 
small baryonic contribution) and Q\ = 0.7. The original 
spectrum of annihilation products is redshifted through 
E' = (1 + z)E, yielding dN(E')/dE' . For gamma rays 
above a few TeV, pair production with the diffuse extra- 
galactic infrared and optical photon backgrounds is an 
important loss mechanism. With neutrinos, such atten- 
uation considerations are unnecessary. (Attenuation in 
Earth would be important for neutrinos above about 100 
TeV, e.g., Ref. [5J], mostly beyond our range of interest.) 

The most relevant factor is the cosmic source intensity 
multiplier, f(z), that reflects the formation history of 
halos. This is a decreasing function of increasing redshift, 
and although the redshift dependence is fairly universal, 
the normalization varies quite considerably, by more than 
an order of magnitude for different profiles [33]. The 
impact of its evolution on the predicted cosmic signal 
is milder once the overall normalization is determined. 
Thus we parametrize it with a simple fitting form 

f(z) = f x io°-9(<-p[-°^]-D-°-i6* } ( 7 ) 



TABLE I: The parameters of Eq. ([T]) for the three profiles 
considered. Scale radius, r a , is in [kpc] and the normalization 
at solar circle, p(R S c), is in units of [GeV cm -3 ]. Estimates of 
the various enhancement factors, J , for the halo (as defined 
in the text) and the normalization of the cosmic intensity 
multiplier, fo, in units of [10 5 ] are also summarized. The last 
row lists canonical values that we use to derive our constraints 
on the DM total annihilation cross section. 
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where /o sets the overall normalization. This normaliza- 
tion is determined by the underlying halo profile, con- 
centration parameter, and inclusion/exclusion of either 
smaller halos or subhalos, and its dependence on the halo 
profile alone can be seen in Table I when the other pa- 
rameters are fixed. The uncertainties of the halo signal 
(including the inner cusp profile, whether a spike or core 
has formed, and substructure), are important for the cos- 
mic signal predictions as well. The cosmological cluster- 
ing of halos, halo mass function and lower mass cutoff are 
additional uncertainties associated with the cosmic sig- 
nal. Thus, the halo signal is overall less uncertain, only 
factors of a few, compared to the cosmic signal, which 
varies by more than an order of magnitude for different 
profiles; compare the entries in Table I. 

The ratio of the DM annihilation halo signal to the 
cosmic signal can be approximated as 

$ C ~ cH^lplf, Jo ' ( ' 

which is independent of (<jav) and m x . Now, using the 
Halo Isotropic signal {Ji so — 0.5) for illustration, this ra- 
tio will be close to unity for /o~5x 10 4 (coincidentally, 
equal to the NFW value as reported in Table I). There- 
fore, even the Halo Isotropic enhancement is significant. 
Its contribution to the apparently diffuse signal observed 
at Earth is comparable to that of the truly cosmic signal 
for profiles as flat as NFW, and it exceeds the cosmic 
signal for profiles flatter than the NFW. This was also 
pointed out in Ref. [34(, specifically for the NFW profile 
with substructures. 

The flatter profiles are well-motivated both theoreti- 
cally and observationally, and hence the Halo Isotropic 
can be an important component. While simulations tend 
to show the presence of a mild cusp, there are also obser- 
vations indicating otherwise, motivating a flatter, cored 
profile p53|. It is not known whether a core is a gen- 
eral feature of the halo structure. In contrast to flatter 
profiles, more cuspy profiles like the Moore profile yield 
much larger fluxes for both the Halo Average and the 
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cosmic signal, and then the Halo Isotropic will be rela- 
tively insignificant. In our calculations below, we use a 
fairly large value of /o — 10 5 to be conservative in com- 
paring constraints that can be obtained from the halo or 
cosmic annihilations. It is unlikely that /o can be any 
larger, unless contributions from the halo substructure 
is significant, which will boost the halo signal as well, 
not significantly modifying their relative strength. Con- 
straints from the halo will be even more important for 
smaller /o values. 

While we use neutrinos for constraining the total anni- 
hilation cross section, the gamma-ray flux has been used 
in many past studies. In this regard, a recent study by 
Ando [33| , which compared the DM annihilation gamma- 
ray signal from the GC to the cosmic signal is of interest. 
It concluded that under the assumptions that annihila- 
tions account for the gamma ray flux from the GC and 
that halo profiles are universal, then annihilations can- 
not be the dominant component of the isotropic extra- 
galactic gamma ray background. That study dealt with 
halos at least as steep as NFW; the Halo Isotropic, on 
the other hand, is most interesting for flatter profiles, as 
we have shown. The Halo Isotropic would increase the 
'apparently' (the cosmic signal plus the isotropic halo 
component) extragalactic gamma ray background, with- 
out violating the constraints placed using the GC data. 



III. NEUTRINO CONSTRAINTS 

To obtain the most robust limits on the total annihi- 
lation cross section for dark matter, we assume that an- 
nihilations produce only neutrinos (the least detectable 
known particles) and compare the resultant neutrino 
signal to the atmospheric neutrino measurements. We 
adopted a simple model in which dark matter particles 
annihilate into pairs of neutrinos, x + X ~~ * v + v, and the 
spectrum of the neutrinos per flavor is a monochromatic 
line with dN v /dE = ^S(E — m x ), where m x is the mass 
of the dark matter particle. The prefactor 2/3 arises un- 
der the assumption that 2 neutrinos are produced per 
annihilated DM pair and that all neutrino flavors are 
equally populated (either in production or through neu- 
trino mixing). This delta function is regularized by the 
redshift integration for the cosmic signal. For the sake 
of a simple graphical representation in Fig. 3, we used 
a narrow Gaussian (width e = m x /20) to represent the 
line signal. 

Strictly speaking, it is not possible to have only neutri- 
nos in the final state, due to electroweak radiative correc- 
tions that lead to Z and W bremsstrahlung (5f|. While 
these corrections do not significantly affect our neutrino 
flux calculations, the vector-boson decay chains do lead 
to gamma-ray fluxes. In turn, these provide a new way 
to model-independently limit the DM total annihilation 
cross section with results comparable to those obtained 
with the cosmic neutrino fluxes [56| . 

In Fig. [3J we illustrate how we set our conservative 
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FIG. 3: The atmospheric neutrino spectrum (u^ + i> M ) and 
the expected increase in the total received intensity from the 
annihilations of a 100 GeV DM particle. While the halo signal 
is sharply peaked in a narrow band of energy, the cosmic signal 
spreads over a wider energy range. To define upper limits on 
the cross section, in each case we use a (different, see Fig. 4) 
cross section such that the annihilation signal would double 
the total received intensity in the displayed energy ranges. 



bound on the total annihilation cross section using neu- 
trinos. The atmospheric neutrino spectrum (specifically, 
v n + ^ju) is derived from measurements by the Frejus, 
Super-Kamiokande and AMANDA detectors, for which 
data are available up to 10 5 GeV 57]. The primary 
source is the decay of pions and kaons produced by cos- 
mic ray collisions with the atmosphere. While the in- 
tensities depend on zenith angle, we adopted an average 
value which will suffice for our purposes (see Ref. [46| for 
details). In the full energy range considered, and in fact 
to higher energies, all of the data are reported as being 
consistent with the theoretically expected flux and en- 
ergy spectrum shape for atmospheric neutrinos. There is 
no evidence of any other signals, e.g., from astrophysical 
sources. 

The DM annihilation signals are calculated using a 
fiducial DM mass of m x = 100 GeV (/ ~ 10 5 for the 
cosmic signal and Ji so — 0.5 for the Halo Isotropic sig- 
nal) and are superimposed on the atmospheric neutrino 
spectrum (per flavor, using + v^). Since the halo sig- 
nal is sharply peaked, we chose an energy bin width of 
Alog 10 i? = 0.3 around E = m x to compare it to the 
atmospheric neutrinos, require it to double the total re- 
ceived intensity in this energy bin, and adjust the cross 
section accordingly. The additional smearing due to red- 
shifting necessitates a larger bin width of A log 10 E = 0.5 
below E — m x for the cosmic signal. These choices are 
within the energy resolution limits of the neutrino detec- 
tors 0,(13. 
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The cross section limits obtained from the above anal- 
ysis are presented in Fig. 2] as a function of the DM mass. 
We show results derived using three different components 
of the halo: the Halo Angular, the Halo Average and 
the Halo Isotropic (shaded regions in order from lighter 
to darker, with each also excluding all larger cross sec- 
tion values), and also the cosmic diffuse signal (dotted 
line). The enhancement factors used in each case are 
listed in Table I. Any larger cross section would pro- 
duce too much intensity and cannot be accommodated 
by the data. We also show the unitarity bound (dot- 
dashed line) [H, [H, [59[ , and the range of cross sections 
for which DM annihilation may flatten cusps (dashed 
line) [H, |42|; cross sections above the lines are excluded. 
The natural scale of the DM annihilation cross section 
(hatched region at the bottom) is shown for a thermal 
relic. However, larger cross sections (not exceeding the 
present limits) are possible, e.g., Ref. [4T| . if dark mat- 
ter is produced non-thermally or acquires mass only in 
the late universe. (For such nonstandard dark matter 
scenarios, it is also possible for the DM to have a large 
scattering cross section with bar yon s, though this is now 
strongly constrained, e.g., Refs. [6fj|.) 

The Halo Isotropic component, which appears in all di- 
rections, is not distinguishable (apart from its spectrum) 
from the cosmic signal, and provides the most robust 
constraint. We also consider the Halo Average and Halo 
Angular components that are relatively more directional, 
yet they can improve on the cosmic constraint by 1-2 or- 
ders of magnitude due to their larger signal intensities. 
If the inner density profile were assumed to be known, 
a greater enhancement in the neutrino signal could be 
possible by restricting to a smaller angular region, as in 
Ref. [ill]. The actual neutrino experiments have differ- 
ent exposures to different parts of the halo, which should 
be taken into account in a more thorough analysis. Thus 
our actual constraint that is indicated by the data is most 
faithfully represented by the Halo Average. 

Several advantages of the halo signal over the cosmic 
signal can be summarized as follows: 

• While both suffer from uncertainties such as the 
concentration parameter and the shape of the halo, 
the halo signal on large scales is overall more cer- 
tain, as noted, and shown in Table I. 

• The isotropic component of the halo signal is es- 
pecially important for flatter profiles, for which it 
dominates over any truly cosmic signal. For cuspy 
profiles, the Halo Angular would be even more con- 
straining than displayed in our Fig. 4. 

• The cosmic signal is broadened in energy by red- 
shifting, making it harder to identify over the 
smoothly varying atmospheric neutrino spectrum. 

• Gamma rays from cosmic DM annihilations are at- 
tenuated at high energies, thus the statement that 
"anything other than neutrinos will be more de- 
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FIG. 4: Constraints on the DM total self-annihilation cross 
section from various components of the Milky Way halo 
(shaded regions excluded). The constraints from the cos- 
mic signal (dotted line) can be compared to that of the Halo 
Isotropic (dark shaded) component. The successive improve- 
ments of the Halo Average (medium shaded) and Halo An- 
gular (light shaded) components can improve on the cosmic 
constraint by 1-2 orders of magnitude. The unitarity bound 
(dot-dashed line) and the cross sections for which annihila- 
tions flatten the cusps of DM halos (dashed-line) are also 
shown; see the text. 

tectable" may not be always fully applicable for the 
cosmic signal (the halo signal will still be present). 

The neutrino limits derived from the halo signal in this 
study not only improve on the previous limit derived from 
the cosmic signal, but also extend their applicability. 

IV. CONCLUSIONS 

Can neutrinos truly be the only final state? Regardless 
of the answer, they can yield the most conservative limit 
on the DM total self-annihilation cross section [46] . Neu- 
trinos are the least detectable particles of the Standard 
Model and all other final states produce more visible sig- 
natures. Thus a limit on their flux (assuming DM only 
annihilates into neutrinos) conservatively bounds the to- 
tal annihilation cross section of dark matter. We focus 
on the importance of the DM annihilation signal aris- 
ing from the halo, while BBM used the truly cosmic 
DM annihilation signal (for decaying DM scenarios see 
Ref. [6l| and references therein). We compare the ex- 
pected spectra of the neutrinos from annihilations to the 
well-measured atmospheric neutrino spectrum and arrive 
at our constraints on the DM self annihilation cross- 
section. We show that, by using the halo signal, the 
previous bounds from the cosmic signal are improved on 
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by 1-2 orders of magnitude. 

This illustrates the point that our intermediate results 
are more generally applicable than just for neutrinos. 
The halo signals have only minor model dependence com- 
pared to the GC and cosmic signals, and also they can be 
significantly larger than the cosmic signals, particularly 
for flatter profiles. These results will also be useful for 
studies with gamma rays. 

An additional enhancement due to possible sub- 
structure of the dark matter or dedicated analyses of the 
data [46| may bring this improvement up to 2-4 orders of 
magnitude. This may turn neutrinos into surprisingly ef- 
fective tools for exploring the full parameter space of the 
total annihilation cross section and mass of the DM, per- 
haps even testing the cross-section scale of thermal relics. 
We emphasize that in the most interesting mass range 
(for lower masses, see Ref. 62]), the cosmic and halo 
neutrino constraints are already significantly improving 
on the unitarity bound [58l 153 ] . Our bounds are espe- 
cially useful at energies > 100 GcV in which there are 
no gamma-ray data on large angular scales, whereas the 
neutrino data go to > 100 TeV. 

In the short term, we strongly encourage dedicated 
analyses of these halo signals by the Super-Kamiokande 
and AMANDA/IceCube collaborations. Our require- 
ment that the DM neutrino signal be 100% as large as 
the atmospheric neutrino background in specified and 
relatively large energy and angular ranges is too conser- 



vative. The atmospheric neutrino background has been 
measured with high statistics over six orders of magni- 
tude in energy, and generally, only much smaller per- 
turbations would be allowed. Together with the other 
strategies for improving sensitivity noted in Ref. [46| , we 
expect that dedicated analyses will lead to much more 
stringent limits than estimated here. 

In the coming age of near-term and proposed super- 
sized detectors, like GLAST, IceCube, LHC and Hyper- 
Kamiokande, together with low background underground 
facilities with unprecedented precision, the days of DM as 
an elusive and hidden substance may finally be coming to 
a conclusion. Neutrinos may play a key role in unveiling 
the true identity of the dark matter. We have demon- 
strated the excellent prospects of neutrino detectors as 
the largest but most conservative tools to constrain the 
DM total annihilation cross section and mass. 
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